A = (1.447391461651648, 2.974068386795948)
B = (0.618583853658536, 0.29905536585365794)
c = Segment(A, B)
C = (3.773956651012925, 0.25368356902715644)
a = Segment(B, C)
b = Segment(C, A)
u: 1.490000015
P: (0, 0, u)
z_1: 1 + (0 * ί)
z_2: ℯ^(((2 * π / 3) * ί))
z_3: ℯ^(((-2 * π / 3) * ί))
f = Ray(P, z_1)
g = Ray(P, z_2)
h = Ray(P, z_3)
α = Angle(z_1, P, z_2)
nr_1, nr_2, F, G, pr_4, I, J, pr_5 = Root(c^(8) + (c^(6) * ((-4 * a^(2)) + (4 * b^(2)) - ((8 * cos(α)^(2)) * b^(2)) - (8 * x^(2)) + ((16 * cos(α)^(2)) * x^(2)) + ((8 * cos(α)^(3)) * x^(2)) - ((16 * cos(α)^(4)) * x^(2)))) + (c^(4) * ((6 * a^(4)) - ((12 * a^(2)) * b^(2)) + (((16 * a^(2)) * cos(α)^(2)) * b^(2)) + (6 * b^(4)) - ((16 * cos(α)^(2)) * b^(4)) + ((16 * cos(α)^(4)) * b^(4)) + ((24 * a^(2)) * x^(2)) - (((48 * a^(2)) * cos(α)^(2)) * x^(2)) - (((8 * a^(2)) * cos(α)^(3)) * x^(2)) + (((32 * a^(2)) * cos(α)^(4)) * x^(2)) - ((24 * b^(2)) * x^(2)) + (((80 * cos(α)^(2)) * b^(2)) * x^(2)) - (((40 * cos(α)^(3)) * b^(2)) * x^(2)) - (((48 * cos(α)^(4)) * b^(2)) * x^(2)) + (((32 * cos(α)^(5)) * b^(2)) * x^(2)) + (24 * x^(4)) - ((88 * cos(α)^(2)) * x^(4)) + ((16 * cos(α)^(3)) * x^(4)) + ((112 * cos(α)^(4)) * x^(4)) - ((64 * cos(α)^(5)) * x^(4)))) + (c^(2) * ((-4 * a^(6)) + ((12 * a^(4)) * b^(2)) - (((8 * a^(4)) * cos(α)^(2)) * b^(2)) - ((12 * a^(2)) * b^(4)) + (((16 * a^(2)) * cos(α)^(2)) * b^(4)) + (4 * b^(6)) - ((8 * cos(α)^(2)) * b^(6)) - ((24 * a^(4)) * x^(2)) + (((48 * a^(4)) * cos(α)^(2)) * x^(2)) - (((8 * a^(4)) * cos(α)^(3)) * x^(2)) - (((16 * a^(4)) * cos(α)^(4)) * x^(2)) + (((48 * a^(2)) * b^(2)) * x^(2)) - ((((128 * a^(2)) * cos(α)^(2)) * b^(2)) * x^(2)) + ((((48 * a^(2)) * cos(α)^(3)) * b^(2)) * x^(2)) + ((((32 * a^(2)) * cos(α)^(5)) * b^(2)) * x^(2)) - ((24 * b^(4)) * x^(2)) + (((80 * cos(α)^(2)) * b^(4)) * x^(2)) - (((40 * cos(α)^(3)) * b^(4)) * x^(2)) - (((48 * cos(α)^(4)) * b^(4)) * x^(2)) + (((32 * cos(α)^(5)) * b^(4)) * x^(2)) - ((48 * a^(2)) * x^(4)) + (((176 * a^(2)) * cos(α)^(2)) * x^(4)) - (((32 * a^(2)) * cos(α)^(3)) * x^(4)) - (((160 * a^(2)) * cos(α)^(4)) * x^(4)) + (((64 * a^(2)) * cos(α)^(6)) * x^(4)) + ((48 * b^(2)) * x^(4)) - (((208 * cos(α)^(2)) * b^(2)) * x^(4)) + (((160 * cos(α)^(3)) * b^(2)) * x^(4)) + (((128 * cos(α)^(4)) * b^(2)) * x^(4)) - (((192 * cos(α)^(5)) * b^(2)) * x^(4)) + (((64 * cos(α)^(6)) * b^(2)) * x^(4)) - (32 * x^(6)) + ((160 * cos(α)^(2)) * x^(6)) - ((96 * cos(α)^(3)) * x^(6)) - ((192 * cos(α)^(4)) * x^(6)) + ((224 * cos(α)^(5)) * x^(6)) - ((64 * cos(α)^(6)) * x^(6)))) - ((-a^(8)) + ((4 * a^(6)) * b^(2)) - ((6 * a^(4)) * b^(4)) + ((4 * a^(2)) * b^(6)) - b^(8) - ((8 * a^(6)) * x^(2)) + (((16 * a^(6)) * cos(α)^(2)) * x^(2)) - (((8 * a^(6)) * cos(α)^(3)) * x^(2)) + (((24 * a^(4)) * b^(2)) * x^(2)) - ((((48 * a^(4)) * cos(α)^(2)) * b^(2)) * x^(2)) + ((((8 * a^(4)) * cos(α)^(3)) * b^(2)) * x^(2)) + ((((16 * a^(4)) * cos(α)^(4)) * b^(2)) * x^(2)) - (((24 * a^(2)) * b^(4)) * x^(2)) + ((((48 * a^(2)) * cos(α)^(2)) * b^(4)) * x^(2)) + ((((8 * a^(2)) * cos(α)^(3)) * b^(4)) * x^(2)) - ((((32 * a^(2)) * cos(α)^(4)) * b^(4)) * x^(2)) + ((8 * b^(6)) * x^(2)) - (((16 * cos(α)^(2)) * b^(6)) * x^(2)) - (((8 * cos(α)^(3)) * b^(6)) * x^(2)) + (((16 * cos(α)^(4)) * b^(6)) * x^(2)) - ((24 * a^(4)) * x^(4)) + (((88 * a^(4)) * cos(α)^(2)) * x^(4)) - (((16 * a^(4)) * cos(α)^(3)) * x^(4)) - (((112 * a^(4)) * cos(α)^(4)) * x^(4)) + (((64 * a^(4)) * cos(α)^(5)) * x^(4)) + (((48 * a^(2)) * b^(2)) * x^(4)) - ((((176 * a^(2)) * cos(α)^(2)) * b^(2)) * x^(4)) + ((((32 * a^(2)) * cos(α)^(3)) * b^(2)) * x^(4)) + ((((160 * a^(2)) * cos(α)^(4)) * b^(2)) * x^(4)) - ((((64 * a^(2)) * cos(α)^(6)) * b^(2)) * x^(4)) - ((24 * b^(4)) * x^(4)) + (((88 * cos(α)^(2)) * b^(4)) * x^(4)) - (((16 * cos(α)^(3)) * b^(4)) * x^(4)) - (((112 * cos(α)^(4)) * b^(4)) * x^(4)) + (((64 * cos(α)^(5)) * b^(4)) * x^(4)) - ((32 * a^(2)) * x^(6)) + (((160 * a^(2)) * cos(α)^(2)) * x^(6)) - (((32 * a^(2)) * cos(α)^(3)) * x^(6)) - (((256 * a^(2)) * cos(α)^(4)) * x^(6)) + (((32 * a^(2)) * cos(α)^(5)) * x^(6)) + (((256 * a^(2)) * cos(α)^(6)) * x^(6)) - (((128 * a^(2)) * cos(α)^(7)) * x^(6)) + ((32 * b^(2)) * x^(6)) - (((160 * cos(α)^(2)) * b^(2)) * x^(6)) + (((96 * cos(α)^(3)) * b^(2)) * x^(6)) + (((192 * cos(α)^(4)) * b^(2)) * x^(6)) - (((224 * cos(α)^(5)) * b^(2)) * x^(6)) + (((64 * cos(α)^(6)) * b^(2)) * x^(6)) - (16 * x^(8)) + ((96 * cos(α)^(2)) * x^(8)) - ((64 * cos(α)^(3)) * x^(8)) - ((144 * cos(α)^(4)) * x^(8)) + ((192 * cos(α)^(5)) * x^(8)) - ((64 * cos(α)^(6)) * x^(8))))
A_1 = Intersect(f, Sphere[P, x(If[u > 1.49, F, If[u < 0.42, G, If[u < 0.71, I, J]]])])
B_1, B_1' = Intersect(g, Sphere[A_1, c])
C_1, D = Intersect(h, Sphere[B_1, a])
p = Plane(A_1, B_1, C_1)
X_{13} = TriangleCenter(A, B, C, 13)
检验数: Distance[C_1, A_1] - b
O = (0, 0)
M: ((x(Intersect[Line[B_1, A_1], Line[Intersect[p, zAxis], Line[A_1, C_1]]]) - x(A_1)) / (x(B_1) - x(A_1)) * (B - A)) + ((x(Intersect[Line[A_1, C_1], Line[Intersect[p, zAxis], Line[B_1, A_1]]]) - x(A_1)) / (x(C_1) - x(A_1)) * (C - A)) + A